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In this paper, we study the effects of parity violation on non-gaussianities of primordial gravi- 
tational waves in the framework of Hofava-Lifshitz theory of gravity, in which high-order spatial 
derivative operators, including the ones violating parity, generically appear. By calculating the three- 
point correlation function, we find that the leading-order contributions to the non-gaussianities come 
from the usual second-order derivative operators, which produce the same bispectrum as that found 
in general relativity. The contributions from high-order spatial n-th derivative operators are always 
suppressed by a factor {H / Mt)"~^ (n > 3), where H denotes the inflationary energy and M« the 
suppression mass scale of the high-order spatial derivative operators of the theory. Therefore, the 
next leading-order contributions come from the three-dimensional gravitational Chern-Simons term. 
With some reasonable arguments, it is shown that this 3-dimensional operator is the only one that 
violates the parity and in the meantime has non- vanishing contributions to non-gaussianities. 



I. INTRODUCTION 

Primordial gravitational waves (PGWs), which are ex- 
pected to be generated during inflation, have attracted 
a great deal of attention recently, as their detections 
would be the direct evidence of inflation, and more im- 
portant the existence of gravitational waves in the uni- 
verse. From the properties of PGWs, such as their power 
spectra and non-gaussianities, we can extract useful in- 
formation about the theory of inflation and gravity. In 
particular, the PGWs produce not only the temperature 
anisotropy, but also a distinguishable signature in the po- 
larization of the cosmic microwave background (CMB) 
[1]. Decomposing the polarization into two modes: one 
is curl-free, the E-mode, and the other is divergence-free, 
the B-mode, one finds that the B-mode pattern cannot 
be produced by density fluctuations. Thus, its detec- 
tion would provide a unique signature for the existence 
of PGWs .2]. 

In addition, PGWs normally produce the TT, EE, BB 
and TE spectra of CMB, but the spectra of TB and EB 
vanish when the parity of the PGWs is conserved Ij. 
However, if the theory is chiral, the power spectra of 
right-hand and left-hand PGWs can have different am- 
plitudes, and then induce non- vanishing TB and EB cor- 
relation in large scales [3]. This provides the opportu- 
nity to directly detect the chiral asymmetry of the the- 
ory by observations [3H5] . Recently, in P| [7] the above 
mentioned problem was addressed in the framework of 
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Hofava-Lifshitz (HL) theory of gravity [S], in which the 
Lorentz symmetry is broken in the ultraviolet (UV) , and 
parity-violating operators generically appear. In partic- 
ular, it was shown that, because of the parity violation 
and non-adiabatic evolution of the modes, a large po- 
larization of PGWs can be produced, and could be well 
within the range of detection of the forthcoming CMB 
observations 0- 

The effects of the parity violation on non-gaussianities 
of PGWs were also studied [HI HO] in the theory with 
the general covariance, and shown that, because of the 
symmetry of the pure de Sitter background, the parity 
violation from Weyl cubic terms have no contributions to 
the non-gussianities, although this is no longer true when 
the coupling of Weyl cubic terms is time-dependent |11) . 
It should be noted that in all these studies the symmetry 
of the general diffeomorphisms of the underlaid theories 
plays an crucial role. On the other hand, in the HL the- 
ory the symmetry is reduced to the foliation-preserving 
diffeomorphisms [8], and the parity- violating operators 
allowed by such a symmetry are quite different from those 
with the general diffeomorphisms. Thus, it is expected 
that in the HL theory some distinguishable features of 
non-gaussianities of PGWs due to these parity-violating 
operators should exist, which may provide a smoking gun 
for the tests of the HL theory in the forthcoming CMB 
observations. With these motivations, in this paper we 
study the non-gaussianities of PGWs in the HL theory, 
and focus ourselves mainly on the effects of the high-order 
spatial operators on the non-gaussianities of PGWs, es- 
pecially on the ones that violate the parity. 

The rest of the paper is organized as follows: In Sec. 
II we first give a very brief review on the HL theory, and 
then restrict ourselves to the model recently proposed in 
[121 113) . where an extra U(l) symmetry is enforced in 
the nonprojectable case, in order to eliminate the spin-0 
gravitons usually appearing in the HL theory. In this sec- 
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tion, we also present the linearized equation of motion of 
the tensor perturbations, originally derived in [7j. In Sec. 
Ill, from the cubic action of tensor perturbations, we cal- 
culate the three-point correlation function and obtain the 
bispectrum of the PGWs, while in Sec IV, we plot the 
shapes of bispectrum produced by both the second-order 
derivative operators and the three-dimensional parity- 
violating Chern-Simons one. In Sec V, we summarize 
our main results. There are also two appendices, A and 
B, in which the cubic action is given explicitly. 

Before processing further, we would like to note that, 
although in this paper we restrict ourselves only to the 
model of the HL theory proposed recently in [T^l [13] , our 
results can be easily generalized to other models, as the 
tensor perturbations are quite similar in all of these mod- 
els [Illll]. In addition, non-gaussianities of PGWs in the 
framework of the general covariant theory with the pro- 
jectability condition was also studied in fl5^, and several 
remarkable features were found. In particular, it was 
found that the terms RijR-i and {V'W^) {^^Rjk) ex- 
hibit a peak at the squeezed limit, while the one RjRl,R^ 
favors the equilateral shape when spins of the three ten- 
sor fields are the same, but peaks in between the equi- 
lateral and squeezed limits when spins are mixed, where 
Rij denotes the 3-dimensional Ricci tensor made of the 
3-dimensional metric gij of the leaves t = Constant, and 
Vi denotes the covariant derivative with respect to gij. 
The consistency with the recently-released Planck ob- 
servations [TB] was also discussed. However, in [TS] the 
parity-violating operators were excluded. Therefore, in 
this paper we shall focus mainly on the effects of these op- 
erators on non-gaussianities, as mentioned above. More- 
over, non-gaussianities of scalar perturbations were also 
studied in the framework of the HL theory, one in the 
curvaton scenario [I7j and the other in inflationary model 
[TSl. and some remarkable features were obtained. 



tioned above can be avoided by properly choosing the 
coupling constants appearing in the theory. Since in all 
of those models, the tensor perturbations are almost the 
same [71 [TH [TSl, without loss of the generality, in this 
paper we shall work with the model proposed in [T^ [T^] . 



A. Action of the Nonprojectable General 
Covariant HL Gravity 



The fundamental variables in the nonprojectable gen- 
eral covariant HL gravity proposed in [T^ [T3] are 



(TV, N\g,,,A,ip), 



where N and denote, respectively, the lapse function 
and shift vector in the Arnowitt-Deser-Misner (ADM) 
decompositions and A and ip are, respectively, the 
U{1) gauge field and Newtonian prepotential (22]. Then, 
the corresponding total action can be cast in the form. 



(2.1) 



II. NONPROJECTABLE GENERAL 
COVARIANT HL GRAVITY AND LINEAR 
TENSOR PERTURBATIONS 



where = l/(167rG) with G being the Newtonian con- 
stant, Lm describes matter fields, and 



By construction, the HL theory is power-counting 
renormalizable [8] . This is achieved by breaking the sym- 
metry of the general covariance in the UV, and include 
only high-order spatial derivative operators, so that it 
remains also unitary, a problem that has been facing for 
a long time in the quantization of gravity [T9|. In the 
low energy, low dimensional operators take over, and it 
is expected that the Lorentz symmetry is "accidentally" 
restored [50]. Since Hofava first proposed it in 2009, the 
theory has attracted a lot of attention, partially because 
of various remarkable features when applied to cosmol- 
ogy [21] , and partially because of some challenging ques- 
tions, such as ghosts, instability and strong coupling. To 
overcome these questions, various models have been pro- 
posed [20], including the ones with an additional local 
U{\) svmmetrv|121 [T51 122 ) . in which the problems men- 
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Here R denotes the Ricci scalar, and Cij the Cotton ten- 



sor, defined by 



Vk{R\-^R5', 



(2.4) 



with e^^^ = 1, etc. A,7„,/3s and Kg are the couphng 
constants of the theory. In terms of Rij, we have |13| . 



]-R^ - ^RRi.R'^ + 3R]RiR'^ + ^RAR 
2 2 
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(2.5) 



where 



G-fc ^ -R^i'S/iR - Ri,\/m'^ - ^R\/^R. (2.6) 

It should be noted that in writing the above action, 
we have excluded all the terms that violate the parity 
[12|, I13j . For our current purpose, we add the fifth and 
third-order spatial derivative operators to the potential 
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ai^3(r) 



aoK,jR,j + a2e'^''RaAjRi) 



(2.7) 



Here the coupling constant ao, Q?i, a2 are dimension- 
less and arbitrary, e^^'^ is the total antisymmetric tensor, 
and W3(r) the 3-dimensional gravitational Chern-Simons 
term ^. denotes the rest of the fifth-order operators 



given in Eq.(2.6) of [13]. Since they have no contribu- 
tions to tensor perturbations, in this paper we shall not 
write them out explicitly. As shown in [7] , bec ause of the 
additional parity violation terms of Eq.(2.7), the non- 



adiabatic evolution of modes lead to a large polarization 
of PGWs, and it could be well within the detection of 
CMB observations, as mentioned above. In this paper, 
we investigate their effects on the non-gaussianities of 
PGWs. 



B. The Linearized Tesnor Perturbations 

The general formulas of the linearized tensor perturba- 
tions were given in [7 , so in the rest of this section we give 
a very brief summary of the main results obtained there, 
in order to to initiate our studies of the non-gaussianities 
of PGWs in the next section. For details, we refer read- 
ers to |7j. Consider a fiat Friedmann- Robertson- Walker 
(FRW) universe. 



N^a{r]), = A = ip^O, 
cjijdx^dx-' — a{ri)'^Sijdx'^dx-' , 



(2., 



where quantities with hats denote the background of the 
FRW universe in the coordinates {ri,x^) — {ri,x,y, z). 
Then, the tensor perturbations are given by, 

6N = 6N' = 6A = Sip = 0, 
6g.ij = a^hijirj^x.). (2.9) 

Assuming that matter fields have no contributions to ten- 
sor perturbations, we find that the quadratic part of the 
total action can be cast in the form. 
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To avoid fine-tuning, a„ and 7„ are expected to be of the 
same order. Then, the field equations for hij read. 



^ To take quantum effects into account, it was proposed to add 
boundary terms AS3 = Y.i l^iM^'^i J^^^^ (fx^O^ into the 
Einstein-Hilbert action at the moment t = tt,, right before the 
inflation started |24| . Clearly, one choice of C is C oc tj(r). We 
thank Jiro Soda for pointing it out to us. 
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where = 1 + "iaoV. / {2M^ a) . 

To study the evolution of hij , we expand it over spatial 
Fourier harmonics, 



(k), (2.12) 



s=R.L 



where Pj^j'(k) are the circular polarization tensors and 



satisfy the relations: ik^e™-' Pif = kp^P^'"''' with = 
1, = -1, and p*"^'^ p^'^"'^ = S'""' [6 . Define uI{tj) ^ 
■^a{ri)Mpiip^{r]) and with the de Sitter background a = 
— 1/ {Hrj) , we obtain 



4(7?) =0, (2.13) 



where 



UJg{k,ri) = a k 1 ~ 6ip^ {eifakr]) + S2{e*ak'i]) 



(2.14) 
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with e* = H/M^ <C 1, and 
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<53^^, ^4^^. (2.15) 
we choose the initial conditions at 
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, =, ul{n,y = ^\h^. (2.16) 

Then, if one assumes that ujs{k, rf) is slowly varying, i.e., 
ujs{k,vi)' 



Q = 



(2.17) 



one can approximatively treat LOs{k,ri) as constant, and 
get the approximate solution of the mode function u^. (77) , 



-i J cjs d7] 



(2.18) 



or, 



2^^il + ^u;,^)e-^f-^^''^. (2.19) 



Promoting i-'f,{ri) to a quantum operator. 



rM^rkas{k) + rJA{~k), 



(2.20) 



one finds that the power spectrum of the tensor pertur- 
bations is given by 



A2 



27r2 



2H^ 



(2.21) 



which has the same expression as that given in general 
relativity. Here as{k) and a|(— fc) are annihilation and 
creation operators, and their commutation relation is 
given by 

[a,{k),al,ik')] = {2Tr)Hss'S{k~k'). 

It should be noted that, in our calculations we have 
assumed (2.17). This condition implies that the adiabatic 
condition is always satisfied (before the modes exit the 
horizon), and thus, there is no important modification 
in the power spectrum of PGWs. Once this condition is 
violated, as shown in [7], some interesting modifications 
on power spectrum and polarization of PGWs become 
possible. For simplification, in this paper we assume that 



(2.17) always holds. 



III. THE INTERACTION HAMILTONIAN AND 
BISPECTRUM 

In this section, we turn to the cubic action and the 
bispectrum of the te nsor p erturbations. The cubic action 
Sg is given by Eq.(A.l|, which can be written in the 
form, 



d-qHiativ)- 



(3.1) 



Then the 3-point correlation function can be computed 
by employing the in- in formalism ,27], 

= -t rdr7'([^^i(?7)V^^^(r7)^^^(r/),Hint(?7')]), 

(3.2) 

where r]i represents the early time when inflation starts, 
and ?7 is a time when the bispectrum is evaluated. A 
good approximation is to extend the integral into the 
whole half axis, 77 S (— cxd,0). After some simple but 
very tedious calculations, it can be shown that the 3- 
point correlation function can be rewritten in the form 



^i{2nf6\ki+k2 + k:i)C' 



TJ/*SlS2S3 



(3.3) 



where Fk^k'k'iv') is given in Appendix B, and 
W^l^llliTj') Is denned as 

xrk:Hri')rkTiv')rkTW)- (3.4) 
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In the de Sitter background, the 3-point correlation func- 
tion reduces to, 



= (2ny6Hk, + k, + fc3)^|^B--3, ^,(3.5) 

Zi n, 2 '''2 "'3 



have 



where 



r>siS2S3 _ \ ^ r P T 



(3.6) 



n=0 



with (5o = 1, and /„ is given by. 



V 

x(l-itj,,77)(l-itj,37/) 1. (3.7) 



Usually, the fc-dependence of the bispectrum receives 
contributions from both the interaction Hamiltonian 
^^int (~ X] ^n) the mode function integration /„. For 
the former, one can see from (3.6) that the high order 



spatial derivative terms do have contributions in bispec- 
trum, but are suppressed by the factor e.^,. Then, the 
leading-order contributions come from the two derivative 
term Fq, which has the same expression as that given in 
general relativity. 

On the other hand, the mode function integration /„ 
in the current case involved very complicated expression, 
and thus it is very hard to perform the integration ex- 
plicitly. In general relativity, to minimize the errors, one 
usually splits the integrals into three different regions: 
one outside the horizon, one around the horizon, and one 
inside the horizon. The contribution from the last re- 
gion vanishes due to the high-frequency oscillate [^H]. A 
similar conclusion is also applicable to the current case, 
because when the high-order derivative terms dominate, 
the oscillation becomes more rapidly than that in gen- 
eral relativity. Thus, using the same arguments, one can 
safely neglect the effects from the high-order derivative 
terms in the last region. 

In the first region, the mode is outsid e the horizon and 
the effective mass term —2/77^ in Eq.(2.13) dominates. 



Then, in this region the corresponding results are also 
the same as those given in general relativity. Therefore, 
the effects from the parity-violating operators come only 
from the region around the horizon. In this region, al- 
though the fc^ term dominates, the high- derivative terms 
still have non-negligible contributions. In order to take 
these effects into account, we expand the integration in 
terms of e* . In particular, to the zeroth-order of e* , we 



Im< i / dri{—ri)" 
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n — 2 i{ki+k2+k3)rj 



X (l-f ikir]){l + ik2v)il + ik3v)j, (3.8) 

which coincides with the mode integration in general rel- 
ativity, as expected. Thus, one immediately obtains the 
bispectrum of the leading-order, 
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K 



K 



(3.9) 



where K = ki + k2 + k^, and which is precisely the bis- 
pectrum of PGWs given in GR. 

Now let us turn to the first order contributions of e*. 
Ignoring all the detailed calculations, it can be shown 
that it takes the form, 



B 



(rvr(^i'^2,fc3) 



Fi + -(Sifci -I- 32^2 -I- S3fc3)Fo 



(3.10) 



Since this term is directly proportional to (5i, from 
Eqs.(2.7) and (2.15) we find that it represents the con- 



tributions of the three-dimensional Chern-Simons term. 

For operators with dimensions n > 4, their contribu- 
tions can be written in the form. 



J, ^ J— 00 



"-4e»^''/„_,_2(?7) 



(3.11) 



where r = (0, n — 2), and fn-r-2{.K, rj) can be expressed 
as 



fn-r-2{K,1]) 



ao + ai{ir]) + 02(177)^ 
+a„_,+i(»?7)"-''+i. 



(3.12) 



here are functions of ki, k2, k^. In particular, the ef- 
fects from fifth-order derivative terms should contribute 
to the bispectrum at the order of e^. But, a careful anal- 
ysis over the integration shows that their contributions 
vanish identically. This result can be easily generalized 
to higher order terms. In fact, for n = 2j + 1 with 
j = 2, 3, 4, the bispectrum of PGWs at the order e"~^ 



always vanishes. This implies that S^py^"^ (fci, /c2, ^3) 
given by Eq.(3.10) represents the only contribution from 



parity violation operators. 

In ref. lOj, Soda et al calculated the bispectrum of 
PGWs from the Weyl cubic terms, and WW^, and 
proved that no contributions from parity violation ap- 
pear in the non-gaussianity of PGWs in pure de Sitter 
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,1.0 




in comparison with the first one in the equilateral limit 
{ki ~ fc2 ~ fca), that is, ~ B+++/81. 

Now we turn to the contributions from the parity- 
violating Chern-Simons term. In this case, because of 
the violation of the parity, for different spin products 
we have four independent configurations. We plot the 
shapes of these four configurations in Figs. ([2]), from 
which it can be seen that all the configurations peak in 

the squeezed limit (Note that for the (+H — ) and ( ) 

cases they peak in the negative direction). More specif- 
ically, we have B(++|(fci, /c2, ^s) = --S(~pv7(^i' ^2, A:3), 
and _Bjtv"(fci,fe,fc3) = --B(pv| 



(fci,fc2,fc3)- 



,1.0 




FIG. 1; Shapes of 

(feifc2fe3)~^B("^"^"^(A:i, fe2, ka 
the equilateral limit. 



{kik2kz) ^Bj+j'py (fci, ^2, fca) and 
). All are normalized to unity in 



background. This is consistent with our current results. 
However, it must be noted that in their considerations, 
the symmetry of the theory is still of the general diffeo- 
morphisms. As a result, only the parity-violating terms 
and WW^ are allowed. These terms are both P-odd 
and T-odd. Thus, when one calculates the bispectrum. 



the two terms produce an integral similar to (3.11 1 with 
n — 2 = 2j — 1 as odd number that is greater than two 
[inilll]- Hence, with the arguments given above, they in- 
deed have no contributions to the bispectrum of PGWs. 



IV. SHAPE OF THE BISPECTRUM 

We are now ready to plot the shapes of the bispec- 
trum. For si = S2 = S3 = 1 and si = S2 = —S3 = 1, 
we plot the shapes of the bispectrum of the leading or- 
der conrtibutions (fcifc2fc3)^^i3^Qpj'^ (fci, /c2, ^3) in Fig. jlj 
Because there is no parity violation in the leading or- 

{ki,k2,k3) and 



der, we have i3^+^+(A;i, ^2, ^3) = B^^^y 
-^(GR) ('^i' ^2, /C3) = B'^ciR){ki,k2,k3). Thus, there are 
only two possible configurations, which both peak at the 
squeezed limit {kr^/ki — 0). As pointed out in [55], 
the second configuration (Fig.( |l(b") |) is sub-dominant. 



V. CONCLUSIONS AND REMARKS 

In this paper, we have investigated the non- 
gaussianities of PGWs generated during the de Sitter 
expansion of the universe in the framework of the HL 
theory, and paied particular attention on the effects 
of the operators that violate the parity. Because of 
the restricted foliation-preserving diffeomorphisms of the 
theory, the parity-violating third and fifth-dimensional 
operators exist generically. By calculating the three- 
point correlation function of the PGWs, we have shown 
that the leading-order contributions still come from the 
second-order spatial derivative terms, and are the same 
as those given in general relativity. The high-order n- 
th spatial derivative terms of the theory also contribute 
to the bispectrum, although their contributions are sup- 
pressed by an factor e"~^. 

More remarkably, we have also found that the three- 
dimensional gravitational Chern-Simons operator u}{T) 
is the only one that violates the parity and meantime 
has non-vanishing contributions to the non-gayssianities 
of PGWs. In comparison with the contributions of 
the second-order operators that produce the same non- 
gaussianity as given in general relativity, its contributions 
are suppressed by the factor e.^,. In addition, operators 
with odd-order and higher than three have no contribu- 
tions to the non-gaussianities of PGWs. 

It should be noted that in obtaining the above results, 
we have assumed that the adiabatic condition (2.17) is 
always satisfied before the modes exit the horizon. If 
this condition fails to hold, the non-adiabatic evolution 
of the modes becomes possible, and hence the integral 
history of the mode function will be dramatically altered. 
Then, large non-gaussianities are expected, although it is 
still an open question how one can extract information 
of PGWs in this case. 
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Appendix A: The Cubic Action of PGWs 

The cubic action for the tensor perturbations can be 
written in the form. 



aM 



a^Mr 
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Then the interaction Hamiltonian is 

-ffint(»7) = J d^xUintiV^x), 



where the even parity terms are 



~ir^2 H Tl~J^i H o 71 JO ^4 



Appendix B: Expression of (r?) 
^feifefesW is given by, 
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n=0 



(A.3) 



F2 



dkihsi + k2S2 + ksssfa'^ {kl + kl + kl) , 

(B.2) 



(B.l) 



Fi = a^dk\kl+4klsi{k2S2 + k3S3)+6kf{k2S2 + k3S3f 



+ Aklsi {kls2 + 4fc2fc^S2 + 4fc|A:3S3 + fcfss) 

+ kf (2kj + Idkjkl + 2kj + 16fc2fc3S2S3 {kl + k^)) 

+ 2klsi (2fcfs2 + Sfca/cssa + 9^2^352 + 9^2^383 + 8^2^352 + 2klsz) 



+ kt 



IQklk^ + 6kl + 2fc3S3S2 (8fc^ + dk^kl + Sfe/cg) 
+ 4fcisi + kjkl + kf) (fc2S2 + ksssf 

+ {ki + kl) (fc2S2 + /e3S3)^|, 

(A:i - 1.2 - h){ln + k2 - k,)(ki - ^2 + k{){ki + A-. + fca) 



(B.3) 
(B.4) 



and the odd parity terms are 



F = 



-adk[kl{Asi - 4s2 - 553) + fc|(4s2 - 4s3 - 5si) + A:|(4s3 - 4si - 652) 
- 2fcjfc^(si - 5s2 - S3 + 2S1S2S3) - 2klkl{si - 4s3 + 2S1S2S3) 



- 2fc2fc3(s2 - 5s3 - si + 2S1S2S3) - 2kik({s2 - 4si + 2S1.S2S3) 

- 2klkl{s3 - 5si - S2 + 2S1S2S3) - 2fc|A;|(s3 - 4s2 + 2S1S2S3) 

+ kfk3{-6si - 5s2 + 6S3 - 4S1S2S3) + kfk2{-3si + 6s2 - 4s3 - 4siS2S3) 

+ k2ki{-6s2 - 5.S3 + 6.S1 - 4S1S2S3) + fc2^3(^3.S2 + 6,S3 - 4.S1 - 4S1S2S3) 

+ A;3fc2(-6s3 — 5si + 6s-2 - 4S1S2S3) + fc3fci(-3s3 + 6si — 4s2 — 4siS2S3) 

+ 2kfklk3{5si + 4s2 - 2s3 + I8S1S2S3) + 2klk2k3{-7si + 7s2 + Tss + 2S1S2S3) 

+ 2fc2fc|fci(5s2 + 4s3 - 2si + I8S1S2S3) + 2A;|fc3/ci(-7s2 + 7s3 + 7si + 2S1S2S3) 

+ 2klklk2{5s3 + 4si - 2s2 + I8S1S2S3) + 2A;|fciA:2(-7s3 + 7si + 7s2 + 2siS2S3)|, (B.5) 



F3 7.5 7.5! I 2 



9i,2r.2jL2 

xi^2kl [k^ - k2klf si + 2kl^kls2 

+8kikl{k2 - ks)kl{k2 + ks)s2 

-klkl {^kls2 + Sk2k3S3 + kl{-5si + 8s2 + 2S3)) 

+klkl(2k3S2 + klkl{-5si + 20s2 - 453) 

+8fc2A;|s3 + 4:klk3S2S3{3si + 4s2 + 253) 

+2fc|(-3si + S2 + IOS3)) 

+kfklkl{8k3Si{si - 83)33 

+ ^252(1 + 8Si(si - S2 + S3))) 

+klklk3Si (sk%S2{-si + S2) 

-8kl + 4fc2fcf S2(4si + 2s2 + 3s3) 

+4fc2^A:2s^(2si +3s2 +4s3) 
+/c|fc| (8siS2 + I6S2S3 - 3) 
+8k2k^ (2 - S2S3 - si(s2 + S3)) 

+fcJ(l + 8S3(-Si + S2 + S3)) 

+2fc|fcf (-1 + 2s3(2(si + S2) + S3))) 

+k'lkl(^4klkj{-si+4s2) 

+2kl{si + 10s2 - 3s3) - 4fc^S3 
+fc|A:|(-4si - 5s2 + 2OS3) 
+A;2A;f S2 (-3 + 8(2si + S2)s3) 
+4A;^fcf (-S2 + 453) 

+2fc|fcf S3(-l + 2S2(2S1 + S2 + 2S3))) 

+klkl A;2fc3(-2si+5s2-8s3) 
+klkl{20si - 4s2 - 5s3) + 2A;f S3 

+ fc^fc|(l + 8S2(S1 + S2 - S3))S3 

+2klk^{10si - + S3) 

+8fc2fc^S2S3(-S2 + S3) 

(-8si - 2s2 + 5s3) 

(6S2 - 9S3 + 8S2S3(-2S2 + S3) 
+ Si(-4 + 8S2S3)) + fc^fcf (-2si + 4S3 + S2(-7 + 8(S2 - 2S3)S3)) 

+ klklk3 k*k3i'iS2 + Si(-7+ 8(si - 2S2)S2) - 2S3) 

+8fc^si + 8fc^^A:^S3 (1 + S2(s2 - S3) - Si(s2 + S3)) 

+8fc2fc3S2 (1 + S3(-S2 + S3) - Si(S2 + S3)) 
+fc|(-4s2 + 6S3 + Si(-9 + 8(S1 + S2 - 2S3)S3)) 

+4fc2fcf (-s3 + 4si) 



+ktkl (2kls3 + klkls3 (1 - 4 + 8si(2s2 + S3)) 
+fc|(-2s2 - 7s3 + 4si(l + 2s3(-2si + S3))) 

+2fc2fc^S2(-l + 2si(si + 2(S2 + S3))) 

+fc^A;|(-9s2 - 4s3 + 2si(3 + 4s2(-2si + S2 + S3)))) 
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